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Abstract 

The dispersionful analogue, by means of Lax formalism, of the zero-genus universal 
Whitham hierarchy together with its algebraic orbit finite-field reductions is considered. 
The theory is illustrated by several significant examples. 

1 Introduction 

In [U [2] Krichever introduced the so-called universal Whitham hierarchies by means of moduli 
spaces of Riemann surfaces of all genera. We are interested in the zero-genus case, which we 
will henceforth refer to as the Whitham hierarchy. This hierarchy contains such particular 
nonlinear integrable systems like dispersionless limits of KdV, Toda, AKNS (or equivalently 
nonlinear Schrodinger) soliton equations. However, the class of integrable systems is much 
wider. There are very important reductions to the Whitham hierarchy, related to the so- 
called algebraic orbits [21 [3], leading to the construction of (l+l)-dimensional integrable 
dispersionless systems that are generated by rational Lax functions. Thus, there arises a 
question of the construction of a theory analogous to the Whitham hierarchy leading to the 
soliton systems, i.e. systems with dispersion. The main reason for the lack of such a complete 
theory are problems with the notion of higher order "finite" poles in the case of algebra of 
pseudo-differential operators. 

The part of the Whitham hierarchy related to the puncture at infinity is nothing but 
dispersionless limit of the well-known infinite-field KP hierarchy. In |1] Krichever considered 
reductions of KP hierarchy generated by rational (pseudo-differential) Lax operators being 
quotients of appropriate purely-differential operators. Besides, he has shown that there exist 
additional symmetries that are not themselves reductions of KP hierarchy. From the more 
general point of view 'quantization' of dispersionless systems from the Whitham hierarchy gen- 
erated by rational Lax functions with "finite" poles of first order is solved in [5]. This solution 
is given by means of Lax operators of KdV type with the additional so-called source terms [6]. 
This class of integrable soliton systems contains, among others, full AKNS hierarchy. In the 
series of articles [7j-[9j Bonora and Xiong investigated integrable differential- difference hier- 
archies of the Toda type arising in the context of the multi-matrix models. They showed how 
this hierarchies, using one flow, can be turned into completely equivalent purely differential 
hierarchies that can be recognized as special rational reductions of the KP hierarchy. They 
also considered dispersionless limit of the above integrable systems. 
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Let us illustrate such a procedure on the simplest example. The Toda hierarchy, being a 
hierarchy of mutually commuting lattice soliton equations, has a Lax formulation expressed 
by shift operators. The pair of linear Lax equations have the form 



n-1 — 



(1.1) 
(1.2) 



where the second equation defines the evolution of the eigenfunction \E'. The compatibility 
between equations (11.11) and (II. 2p leads to the Toda equation. It is known that the above 
auxiliary linear problem can be rewritten with the use of pseudo-differential operators as it 
follows from (Q that = {d^ + w„)^'„ and = {d^ + w„_i)-^^'„ [3, [ID] -[IS]. Then, 
the linear equation (II. ip takes the form 



of the eigenvalue problem related to the AKNS system. Consequently, the full Toda hierarchy 
in the representation of pseudo-differential operators takes the form of the well known AKNS 
hierarchy. Following the above idea in the opposite direction one is able to construct, with 
the aid of shift operators, a dispersionful analogue of the Whitham hierarchy. 

Very recently Takasaki considered dispersionless limit of mult i- component KP hierarchy 
with charges (lattice variables), that can be treated as a generalization of the Toda hierarchy 
|13j . In this sense the original Toda hierarchy is equivalent to the two-component KP hier- 
archy. Rewriting the mult i- component KP hierarchy in the scalar Lax formalism Takasaki 
showed that the above dispersionless limit is the Whitham hierarchy [13]. Further in [H] 
Takasaki and Takebe analysed the above issue from the point of view of Hirota equations and 
Fay identities. 

The theory of Whitham hierarchy and its dispersionful analogue is inseparably connected 
with problems of systematic constructions of integrable nonlinear dynamical systems. It is 
well known that a very powerful tool, called the classical i?-matrix formalism, can be used for 
systematic construction of (l+l)-dimensional field and lattice integrable dispersive systems 
(soliton systems) [13]- [IH] as well as dispersionless integrable field systems yil|-L22j. More- 
over, the i?-matrix approach allows a construction of Hamiltonian structures and conserved 
quantities. 

The aim of this article is to consider the dispersionful analogue of the Whitham hierarchy 
and its algebraic orbit reductions on the level of the explicit Lax representations. The article 
is divided onto two parallel parts concerning dispersionless and dispersive cases, respectively. 
First the basic facts about the Whitham hierarchy together with its algebraic orbits, given by 
some meromorphic functions, are presented. Next an alternative formulation of the original 
Whitham hierarchy is constructed. This formulation is the one that in the formal quantization 
procedure yields the dispersionful Whitham hierarchy. Initially the dispersionful Whitham 
hierarchy is presented by means of pseudo-differential operators, but the one-to-one analogue 
of the original Whitham hierarchy is given with the simultanous use of the shift operators 
following reference [13]. This is a multi-component Toda-like formulation of the dispersionful 
Whitham hierarchy. The finite-field reductions being counterparts of algebraic orbits are also 
considered. The whole theory is illustrated by the rich set of significant examples. Finally 
some conlusions are given. 
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2 Whit ham hierarchy 

2.1 Basic definitions and Hamilton- Jacobi problems 

Consider a set of Lax functions Za{p,t) {a G {oo, 1, A^}) in a complex variable p being 
formal Laurent series at oo and 'finite' punctures Qi on the Riemann sphere 

oo 

Zoo = P+^Uoc,l{t)p~^ (2.1) 



1=1 



oo 



2, = Mi,_i(t)(p-g,(t))-i + ^Mi,z(t)(p-gi(t))' zel,...,iV, (2.2) 

1=0 

respectively. In fl2.1tl2.2p the coefficients Ua,i and the poles qi are smooth dynamical fields 
depending on the infinite set of evolution parameters ('times') t = {tan}- One of these times 
is distinguished as a spatial variable: tool ='■ x. Times tan are coupled with the generating 
functions VLan defined in the following way 

for a = OO, n = l,2,3, ... , 
-(22)|„,_) fora = z, n = l,2,3,... , (2.3) 
In [p — Qi) for a = 2, n = 0, 

where {■){a,+} are the projections onto the principal parts of Laurent expansions at poles 
oo and qi such that (Xlfc akp''){oo,+) ■= Y^k^o cikP^ and (Xlfc CLk{p-qiY){i-) ■= Y.k<o <^k{p-qiY, 
respectively. Then, the Whitham hierarchy is the set of zero-curvature equations 

^-|^ + {(]_,(],„} = (2.4) 

for the following hierarchies of mutually commuting evolution systems given in the Lax form 

dZa 

— = {Vtjjmi Za] ^ (2.5) 

Otpni 

where the Poisson bracket is the canonical one 

{f^d} dp dx dx dp ^ 

Let us use the short notation tam ^ — ^ t/3n for the equation, related to fixed times tam and 
tffm from the Whitham hierarchy (12.41) . Notice, that VL^oi = p thus tan < — > ^ooi is trivially 
satisfied and the equations from Lax hierarchies (12. 5p for = oo,m = 1 are the translation 
symmetries {za)taai — {za)x in the variable x chosen as a spatial one. 

The above theory of Whitham hierarchy is well defined due to the standard approach 
following from the classical i?-matrix formalism [21 [ISl 122] • The Lax hierarchy (12.51) for 
q; = /3 = oo is the well known infinite-field Benney chain, containing many well known finite- 
field reductions being dispersionless limits (see for example [21]) of soliton systems [I6] from 
the KP hierarchy. 

Examing the bi-Hamiltonian structure of equations (12. 5p in the spirit of the classical R- 
matrix theory [22], one observes that hierarchies (12. 5p for different fixed a and n ^ 1 form 



different bi-Hamiltonian hierarchies, mutually commuting, with respect to the same Poisson 
tensors. The symmetries with n = 1 are the starting symmetries for the recurrence procedure. 
In the same way equations for a = i and t?, = 0, constructed for Q = ln(p — g,), generate other 
bi-Hamiltonian hierarchies, which except the first symmetries are not constructed within the 
above scheme of the Whitham hierarchy with generating functions restricted to (12.31) . In 
this sense the above construction is still incomplete. The appearance of the logarithmic 
function in (12.31) could be explained in several ways. For example in [23] all terms (12. 3p . 
including the logarithmic ones, are constructed by the Lie algebraic splitting of Hamiltonian 
vector fields. How to generalize the Whitham hierarchy so that it contains the whole bi- 
Hamiltonian hierarchies generated by functions containing logarithmic terms is not clear yet. 
Restricting to the particular reductions of the Whitham hierarchy given by meromorphic 
Lax functions with only two poles, first oo and second 'finite' one, one can manage with 
this problem proceeding like in pi], where logarithmic terms were represented by a sum of 
two infinite series convergent at these poles. For the same class of Lax functions the more 
correct approach is given in [25], where the logarithmic terms are constructed using contour 
integrals. However, the generalization of such approaches to the whole Whitham hierarchy 
makes problems. 

Example 2.1 Consider the part of the Whitham hierarchy involving times: tool ='■ too2 ='■ 
t, til ='■ Hi one? tio =: Si. The above notation will be used in the whole article. The functions 
(12.31) related to these times have the form 

Qooi=P VLoo2=P^ + '2w VLh = -ai{p - qiY^ fiio = ln(p - g^), 

where w := Uoo,i, Oi := Then, for t < > yi one finds 

{qi)t = {2w + qi)x {ai)t = {2aiqi)x Wy^ = {ai)^ (2.7) 

the generalized Benney gas system. For yi < — > yj and Si < — ^ yj we have 

and 

(2.9) 

fe).,: ^^3. (2.10) 

respectively. Eliminating qi from (12. 9p one obtains the dispersionless multi- dimensional Toda 
equation or equivalently the Boyer-Finley equation, which appears in general relativity theory. 
For Si < — > t one finds 

Ws, = {qi)x 



ili - (lj)x 
Qi - Qj 



and the first equation from (12. 7p . For Si < — > Sj we have the first equation from (I2.10p . The 
set of equations from this example has been considered in 12^ . 



Furthermore, there exists a set of pseudo-potential functions Sa{za, t) satisfying a set of 
the Hamilton- Jacobi problems 

= Za{dxSa,t) (2.11) 

fi^J5,§«,t) (2.12) 



dt 



and compatibility conditions 



/3n 



- (2.13) 



since p = dx^a- The above equations must be correctly understood, i.e. (12. lip are Hamilton- 
Jacobi equations with fixed 'energy' levels and therefore §>a depend on Za treated as pa- 
rameters. Eliminating dxEi from (12.11112.121) the compatibility conditions (12.131) lead to the 
Whitham hierarchy (12. 4p and simultaneously to Lax hierarchies (12.50 . 

2.2 Algebraic orbit reductions 

There is important class of finite-field reductions of the Whitham hierarchy, being the so-called 
algebraic orbits given by meromorphic functions E satisfying 



E = = z^' for Ha e Zh 
Then, the most general form of E is given by 



noo—2 N rii 

E = p^-+J2 «ooy + E E ""^'i^^P - l^y'^ (2.14) 

1=0 i=l li=l 

and then the dynamical fields from functions z^ are given by polynomials of fields from 
(I2.14p . i.e. expanding E at oo E = z^ and expanding E at Qi E = z^\ The goal of such 
reductions of the Whitham hierarchy is the construction of (l+l)-dimensional integrable 
finite-field dispersionless systems. By integrable systems we understand those which have 
infinite hierarchy of commuting symmetries. In this case pseudo-potential functions reduce to 
one Sq(2q,, t) := S(i?, t) and the set of Lax hierarchies (12. 5p reduces also to one Lax hierarchy 
in the form 

dE 

— = {VtjSm-, E} , 

OlfBrn 

where the functions VLpm are generated around poles of E. 

A more general theory of meromorphic Lax representations, not only for the canonical 
Poisson bracket (12. 6p . including the case (I2.14p allowing a construction of integrable disper- 
sionless systems together with multi-Hamiltonian structures is presented in [22] . In the same 
paper, it is also shown that one can construct Lax hierarchies not only at poles of mero- 
morphic Lax function E but also at its zeros. Another class of reductions of the Whitham 
hierarchy related to algebraic orbits is presented in [3J. 

In all examples in this article we will in general present only the first nontrivial evolution 
equations from related hierarchies. 



Example 2.2 The dispersionless AKNS hierarchy. 

We will consider the meromorphic Lax function in the form 

E = p + '^ttiip- Qi)'^ = = Zi, (2.15) 

i 

i.e. the case of fl2.14p for n^o = = 1. Then the system generated by f2oo2 = p"^ + 2 Oj is 

{ai)t = 2{aiqi)^ 

(2.16) 



{q{)t = 2i^a,+qi^ 



This is the dispersionless limit of multi- component AKNS system. For VLko = ln{p — qk) we 
have the following equations 



ai 



{(li)sk = (ln| gi - Qfcl) 



X 



k^k 

{qk)s, = (Inafc)^. , 



where i ^ k. For i^^i = ~ak{p — qu) one finds 

aiOk 



X 

ak 



Qi - Qk 



X 



(2.17) 



i^k 



{qk)y, = {qk)x - ( 

i^k ^ 



qi - qky ^ 



where i ^ k. The equations from this example are of course compatible with the equation from 
Example \2.1\ as in this case w = J2i ^j- 

Example 2.3 Consider Lax function (I2.14p for N = 1 with the pole of first order at oo and 
second order pole at q := qi: 

E = p + u{p- q)-^ + v{p - q^ = z^ = z\. (2.18) 

Let y := yi, s := Si. Then, one finds the following equations 

Ut = 2{uq)x + 2vn, 

VLao2=P^ + 2u =^ vt = 2v^q + Avqx (2-19) 

qt = 2u^ + 2qq^, 



Us = Qx 

UVx 

fiio = ln(p-g) =^ Vs = Ux-— (2.20) 



V. 



X 



Uy = iVv)x 

^n = -Mp-q)-' =^ "^^^^l^^iV; (2.21) 



u 



^V^Jx' 



Now w = u from Example \2.1[ 



Example 2.4 Let us take now (12.141) for N = 1 with rioo = 2 and rii = 1, i.e. 

a 



E = + u-\ = zl^ = zi. 

p-q 

In this case E = floo2 — ^ii, thus we have that Ef = {Qoo2,E} = {Qii,E} = Ey (y := yi) 
and consequently y and t can be identified, i.e. t = y. The system for these times is 

Ut = 2ax 

n^2=p'^ + u =^ at= (2ag)^ (2.22) 

qt = {q^ + u)x. 

For the times r := toos one? s := Si one finds the followinq systems 

Ur = ^^^^ + 3aQ 

3 3 / 3 \ 

fioos = P^ + -up + -a =^ Or = \^3aq'^ + -juaj (2.23) 

3 3 
q^ + -uq + -a 



and 

Us = 2q^ 

VtiQ = \n{p - q) =^ Os = {q^ + u)x 

qs = (Ina)^. 

2.3 Alternative formulation 

The choice of the evolution parameter tool to be considered as a spatial parameter x is 
ambiguous. Since flko = ln(p — qk) by (12.121) the following relation is valid 

= exp {ds,S^) + qk, (2.24) 



where Sk ■= tko- Now fixing k, we can eliminate dsf.E>a from fl2.im2.l2|) instead of 9^Sq and 
consider Sk as a spatial variable. Then, for = exp{dsf.Sa) the Whitham hierarchy (12.41) and 
related Lax hierarchies (12.51) have the same form, but the Poisson bracket is given by 

Sf \ = x f^^_^^\ /2 25) 

^ \d\kdsk dskdXkJ ' 

In such a situation the variables p and should be considered rather only as auxiliary 
parameters. On the level of zero-curvature equations (12. 4p and Lax hierarchies (12. 5p the 
above procedure is little bit more complicated. /^From (I2.24p it follows that at first, one has 
to perform the transformation X^. = p - q^. Hence, dp h-^ c^a^ and dt^„ dt^„ - {qk)t^^dx^. 
Next, all terms in (12.41) and (12. 5p containing derivatives with respect to x have to be replaced 
using equations t^n < — > Sk for VLko = InA^ and (12.51) for /3 = /c, m = 0. As a result, one 
obtains, in a preserved form, the Whitham hierarchy (12. 4p and related Lax hierarchies (12.50 
for the Poisson bracket given by (12.250 . On the level of explicit equations both formulations 
give a compatible set of equations. In the case of explicit equations from Lax hierarchies (12.50 
the above passage between both formulations relies on the change between spatial variable x 
and Sk through use of explicit equations (12. 5p for /3 = A;, m = or = oo, m = 1, respectively. 

Of course in the new auxiliary variable A^ the Lax functions (12.1112.20 take appropriate 
form of Laurent series. The "finite" poles qi are then moved to qi — qk- Particularly Zi takes 
the form of a Laurent series at 

oo 

Zi = fti.-iAj ^ + 0'i,iXi\ 

1=0 

so the dynamical field qi in Zi is missed. Then, the evolution of qi can be calculated from the 
representation of 2:00 as 



Xk + Qk + 0'(x.,iXk ■ 



00 

1=1 



In fact, we can use all auxiliary variables p and Aj (for all i) simultaneously. This formu- 
lation of the Whitham hierarchy is correct if all equations (relations) are understood on the 
level of the Hamiltonian-Jacobi problems, i.e. p = dxE>a, A, = exp{ds3a) are not considered 
as independent variables but as shortened notation. Then, the Poisson bracket is given in the 
following general form 

If j = _ + \^aY— — - — — 

dp dx dx dp ^ * V ^\ 9si dsi dXi 

i 

In this case one has to first evaluate the Poisson bracket and then use the relation p = Xi + qi = 
Xj + qj to draw consistent equations from (12.40 or (12.50 . However, one has to proceed carefully 
in order to not lose the evolution of the fields q^. 

Example 2.5 With the use ofXi, the functions (12. 3p related to the set of evolution parameters 
X, t, Si and i/i are given in the form 



^ooi = Xk + qk fioo2 = Xl + 2qkXk + ql + 2w Qn = -aj(Aj) ^ Qio = In Aj, 



where k is arbitrary from allowed values. Then, one obtains a set of equations being consistent 
with the set of equations from Example \2.1[ 



2.4 Finite-field reductions 

The meromorphic reductions fl2.14p can be represented using variables Xi = p — qi simultane- 
ously in the form 

noo-2 N rii 

E = \l-+ uXl--' + J2 ^ooAk + E E V\ (2.26) 

1=0 i=l k=l 

where k is arbitrary and u = riooqk- However, as we are interested in the construction of 
(H-l)-dimensional equations we should use the form of E only with one auxiliary parameter 
Afc for a fixed k. In this case, Qko = In Xk leads to the trivially satisfied equations related 
to the translational symmetry, this time, in the variable Sk- Therefore, the equation from 
the hierarchy (12.51) for a = A; which was earlier calculated for Qko = — Qk) can now be 
calculated for ^ooi = Xk + Qk but in the reversed form, i.e. with x as an evolution parameter. 

Example 2.6 Multi- component dispersionless Toda hierarchy. 

In the auxiliary parameter Xk = p — qk the Lax function f l2.15l) takes the form 

E = Xk + qk + dkXk + ^ ai{Xk + qk- qi)'^ ■ 

The first nontrivial system is a counterpart of translational symmetry in x for (12.150 

{ai)x = {qi - qk){ai)s^ + ^Mklsk 
^ooi = Xk + qk ^ fe)x = (g.-g.)(?.)., + 5^(a,).,. (2-27) 

This is the multi- component dispersionless Toda equation. The counterpart of dispersionless 
AKNS system (12.161) is only the next system from the hierarchy calculated at oo. For Vtki = 
—UkX^^ one obtains a counterpart of the system (I2.17p . which can be also calculated directly 
from (1^371) with the use of (if^ . 

\ai)yk — - - ^ T ^ 

qi — qk 

{ak)y^ = ak{qk)s, - ^{aj)y^ 

As a result, the dispersionless AKNS hierarchy from Example \2.2\ with x as a spatial variable, 
becomes a dispersionless multi- component Toda hierarchy with the spatial variable Sk- 

Example 2.7 In X := Xi = p — q the Lax operator (\2.18h has the form 

E = X + q + uX-^ + vX-^ (2.28) 

and the first nontrivial equation becomes 

Ux = uqs + Vs 

^ooi = X-\-q =^ Vx = 2vqs (2.29) 

Qx = Us- 



J- \J 



The counterpart of fl2.20p is only the next system from the hierarchy. Then, 



Uy 








Vy 











is the counterpart of fl2.2ip . 



Example 2.8 In terms of X := p — q we can also rewrite the algebraic orbit of Example \2.4 
as: 

E = \^ + 2q\ + + U + aX'^ . 



So, we have 



Ux = as- qus 

n^i = X + q =^ = (2.30) 



As before E = Q002 — ^11 ■ The counterpart of (12.221) is 

Ut = 2aqs 

= X"^ + 2qX + q"^ + u =^ at = a{q'^ + u)s (2.31) 

qt = tts- 

3 Dispersionful Whitham hierarchy 

3.1 Pseudo-differential operators and auxiliary linear problems 

In this section we will use the same set of evolution parameters t^n and the same notation 
X '■= tool, Si := tiQ and yi := tn as in the previous section. 

It is well known that the Benney momentum chain can be obtained from the infinite-field 
KP hierarchy in the quasi-classical limit, which suggests that the dispersive counterpart of 
the Whitham hierarchy can be constructed by means of pseudo-differentials operators d := dx 
satisfying the generalized Leibniz rule 



d"^uix) = J2 uix^d""-^, (3.1) 



where {'^) = (-l)"(""'t""^) < 0. From 1^ it follows that 

((9 — v)u = u{d — v) + Ux 
{d — v)^^u = u{d — v)^^ — (9 — v)^^Ux{d — v)^^ 

= u{d - v)-^ - Ux{d - v)-^ + U2x{d - v)-^ - ... , 



where v can be equal zero. Some other useful formulae, needed in further calculations, are 

d{d-v)-^ = l + v{d-v)-^ 
{d-v)-^d= l + {d-v)-^v 

{d-v)-^u = u{d-v + ^y^ 
[(d-v)-']^ = id~vr\id-v)-\ 

Let us consider the following algebras, generated respectively by elements d, and 

(d-qd,id-q,)-\ 



3oo = 0(oo,+) © 0(oo = < 5Z 

Qi = 0(i,+) © 0(i,-) = \ QiY \ ® l^ai{d 




1^0 

in which the operation, with respect to the above rules, is associative. The Lie algebra struc- 
ture is defined through the commutator [A, B] = AB — BA, where A, B & q^, and these 
algebras are decomposed into Lie subalgebras. We would like to have on the 'quantum' level 
of pseudo-differential operators a theory parallel to the complex calculus of meromorphic func- 
tions, but after quantization we do not have to our disposal such a crucial tool as expansion 
into Taylor series and the clear-cut notion of poles for pseudo-differential operators. However, 
it is possible to construct formal rules of expansions which are similar to the expansions for 
meromorphic functions, i.e. 0oo,0j C 0(j,+) for j i and 0(i,±) C 0(oo,±); where a C b means 
that all elements from a can be expanded as a series from b. Hence, the reasoning is analogous 
to the one on the Whitham hierarchy level. This is a crucial fact for the consistency of the 
equations (13.71) and (13. 8p . The above scheme will be considered in detail in the forthcoming 
article |2g. 

We define the infinite-field Lax operators in the form 

oo 

L^ = d + Y,'^oo,i{^)9-' (3.2) 
1=1 

oo 

L, = ii,,_i(t)(a - qi{t))-' + ^ uu{t){d - q,{t)y i G 1, iV, (3.3) 

/=o 

where all coefficients depend on the same set of evolution parameters as before. The first one 
Loo defines the well-known KP hierarchy. Then 

l].„(t):=|^^-)(--) '°^" = ^ (3.4) 

fora = z, 

where n = 1,2,3,.... Due to the lack of a good counterpart of logarithmic function, formulated 
by means of pseudo-differential operators, we omitted the case related to VLiq in (12. 3p with the 
natural logarithm. It will be considered separately. 

Now we can define a set of auxiliary linear equations 

L„^« = Zo,^^ (3.5) 



dt 



i3n 



(3.6) 



for some eigenfunctions \E'a(t) and eigenvalues Za- Then, the compatibihty conditions between 
(13.51) and (13.61) lead to the Lax hierarchies 

[dt^^-npn,K]^a=^. (3.7) 

The compatibility conditions among equations (13.61) lead to the zero-curvature equations 

[d,_-n^^,dt,^-np^]^, = Q (3.8) 

of Lax hierarchies (13. 7p . They are written as operation on eigenfunctions as this form guaran- 
tees that they are invariant with respect to multi-component Toda-like formulation presented 
in Section ing Still we have to argue that they are well defined. One can show that equations 
(13.71) are self-consistent, i.e. Lax operators (I3.2tl3.3l) are well normalized and the right- and 
left-hand sides of Lax hierarchies can be written in the same form. The related i?-matrices 
are defined through decomposition of Lie algebras Qa and importantly they are mutually 
consistent due to the above rules of expansion. The mutual commutativity between symme- 
tries from fixed Lax hierarchy (13. 7p follows now from the scheme of i?-matrix theory [THi ITS] . 
We will use the notation t^m ^ — ^ tpn for the related equation from the zero- curvature hi- 
erarchy (13. 8p . Notice, that fioo,i = d- Hence, tan ^ — > ^ooi is trivially satisfied and in the 
representation of pseudo-differential operators {La)t^^ = {La)x- 

The so-called quasi-classical limit of linear equations (I3.5ti3.6p leads to the Hamiltonian- 
Jacobi equations (I2.11H2.1"21) for n 7^ 0. To see that, at first one has to perform the transfor- 
mation 

tan ^ jtan =^ dt^„ ^^ Hdt^„, 

where H is treated as a deformation parameter. Next, assuming the WKB form of the eigen- 
functions and that dynamical fields have quasi-classical counterparts, i.e. 

^„ = exp + Oih)^ 

one takes the limit h 0. Then, the fields Ua,i{t) can be identified with those from Lax 
functions (12. 1112. 2p . Hence, the zero-curvature equations (13. Sp can be considered as dispersive 
analogues of the Whitham equations (12. 4p without the case related to n = 0. However, such 
an approach with the use of the formalism of pseudo-differential operators has fundamental 
disadvantages. The first one is that the coefficients of for tt, ^ 2 do not have 'limited' 
or 'compact' form, which means that they are differential polynomials in all, infinitely many, 
coefficients from Lax operators Lj. It makes particular problems when one tries to define 
counterparts of algebraic orbits for the original Whitham hierarchy. 

Example 3.1 Consider equations from the dispersive Whitham hierarchy related to the times: 
tool ='■ X, too2 ='■ t, til ='■ Ui for which the related operators (13.40 have the form 

where w := Moo,i, Oj := Ui^i. Equations for Qooi ore satisfied trivially. Then, for t < > yi 

one finds 



{ai)t = {ai)2x + 2{aiqi)^ {qi)t = -{qi)2x + '^qMx + '^w^ Wy^ = (ui)^, (3.9) 



i.e. the dispersive counterpart of the generalized Benney gas system (12.71) . To calculate the 
equation related to yi < — > yj it is better to rewrite Qn in the form fia^a = —'^id~^(j)i'^a, 
where ai = i/jicpi and qi = — (In^j)^;. Then, the dispersive versions of (12. 8p are 

where means the formal integration operation which cannot be confused with . With 
the use of ifji and (pi the system (13. 9p becomes 

{i)i)t = ii^i)2x + '2w^pi {(j)i)t = -(0i)2x - 2«;0i Wy^ = (^i0i)x. 

3.2 Finite-field reductions 

After quantization we would like to construct reductions of the above dispersive Whitham 
hierarchy being counterparts to algebraic orbits for the Whitham hierarchy, i.e. we are looking 
for a Lax operator £ satisfying the following linear equation 

= * = i^r* = noo, ^ 1, (3.10) 

where now the eigenfunctions \E'q, reduce to one \E' and we have only one eigenvalue E. The 
related Lax hierarchy takes the form 

[dt,„-npn,C]'^ = 0. (3.11) 

The disadvantage mentioned above practically makes it impossible to construct C related to 
Ui ^ 2 with the use of pseudo-differential operators. For = 1 this operator can be defined 
as 

(noo-2 N \ 

+ J2 uid' + J2aiid- qi)~' ^. (3.12) 
1=0 i=l ) 

It is useful to replace ai{d — qi)~^ by the operator '0jc)~"'^0j, where = ^/'i^j and qi = — (In^j)^. 
Then, (13.121) takes the form 

/ noo-2 N \ 

= a"- + ^ uid' + J2 ^id~'(Pi ^. (3.13) 

V 1=0 1=1 / 

The Lax operators in the form (I3.13P have been considered in [6] as reductions of the KP 
hierarchy, i.e. only for a = oo. In a more general context, i.e. for all a, the theory of such 
operators together with Hamiltonian structures is presented in [5]. The fields ipi, <pi are the 
so-called source terms as i/ji and (pi are eigenfunctions and adjoint-eigenfunctions of the Lax 
hierarchy (13. lip , respectively, i.e. ipi satisfy linear equations (13.60 on eigenfunction ^. 

Example 3.2 The multi- component AKNS hierarchy. 

Consider the case = 1 of (13.120 or equivalently (13.130 . that is 

C^= (^^ + J2^^i^-q^r'^'^ = ^a + ^V'^^-V^j ^- (3.14) 



Then for fioo2^ = {d"^ + 2 aj)\l/ = {d"^ + 2 one calculates the multi- component 

AKNS system in a form similar to its dispersionless limit (12.161) 

{ai)t = iai)2x + '^{caq^x 

{Qi)t = -{qi)2x + 2 ^{aj)x + 2qi{qi)^ 
j 

or in the representation with source fields 

j 

For Qki"^ = —'4'kd~^ 4>k^ one finds a dispersive counterpart of f l2.17p 

jf^k 

My, = {4>k)x + ^(f)jd;\ijj(f)k), 

ii^k 

where i ^ k. Obviously the above systems are compatible with the equations from Example \3.1\ 
as w = 



Example 3.3 Let us take dispersionful analogue of the Lax function from Example \2.4\ then 
we have 

/:^= (^d'^ + u + a{d-q)-^)^. 
We have that = (fioo2 — ^u) ^- So, as before Ct = Cy, i.e. t = y. Thus, 

Ut = 2ax 

J^oo2^ = ((9^ + u) ^ =^ at = {ax + 2aq)x 

Qt = {-Qx + q^ + u)x, 

which is the dispersionful version of (I2.22p . For fioos^ = {d^ + |w(9 + + 2a)) \1/ one 
obtains the dispersionful version of fl2.23p 

/^l 3 3 2 o 
Ur=[ |M2x + -ax + + 3aq 

2 3 

a2x + 'icLxQ + ^aq + -ua 

3 3 3 

q2x - "iqqx - -Ux + q^ + -uq + -a 



3.3 Multi-component Toda-like formulation 

In [IH] a unified approach to systematic construction of field and lattice soliton systems was 
presented. It was shown that Lie algebras of pseudo-differential operators can be consid- 
ered as Weyl-Moyal-like deformations of algebras with Poisson structure given by canonical 
Poisson bracket f l2.6p . The Weyl-Moyal-like deformation is a special case of the deformation 
quantization. On the other hand, Lie algebras of shift operators in the same procedure can 
be considered as deformations of algebras with Poisson bracket in the form fl2.25p . which 
suggests that shift operators and similar relation to (11.21) can be helpful in solving problems 
appearing in the pseudo-differential representation of dispersionful analogue of the Whitham 
hierarchy. Hence, instead of constructing the operator counterpart of the logarithm functions 
Qio we will use, following [13], the 'quantum' analogue of the relation f l2.24p . given at this 
moment ad- hoc as 

^ = (e^^. + (3.16) 

where e^^'if = Xli^o is the contracted formula for Taylor expansion. Then, the quasi- 

classical limit of (I3.16P is exactly the relation (I2.24p . Let us denote by := e^^fc the shift 
operators that when acting on dynamical fields satisfy 

A^M(sfe) = u{sk + m)A™ meZ. 

The field qk has shifted the argument in (13.160 . so then it follows that 

Afc^, = {d- Qkisk + 1))^„ =^ A, = {d- qk)-'^a. (3.17) 

/^From (I3.17P we have also the relation between different shift operators 

Ai^^ = (Afc + Qkisk + 1) - Qiisi + 1)) (3.18) 

It is important to remember that the above rules are valid only when acting on eigenfunctions 
"^a- Other relations which may be useful are 

Ak{Ak-v)-^ =v{Ak-v)'^ + l 
{Ak-v)-'Ak = iAk-v)-^v + l 

(Afe - vy^u = u{sk - l)(Afc - v)-^ - (Afc - vy^v [u - u{sk - 1)] (A^ - v)'^ 
(A, - v)-'u = u{sk - 1) (a, - v^^y 
[(A, - v)-']^ = (Afc - vy'vtiAk - v)-\ 

We can add the linear relations (I3.16P to the auxiliary linear equations (13. 6p and complete 
the set of zero-curvature equations (13. 8p by 

Si " W : [d - A, - qi{si + 1), dt^^ - n^rr.] ^7 = (3.19) 
Si< >Sj: [d-A,-q,{s, + l),d-A,-qj{sj + l)]^^ = 0. (3.20) 

Adequately the Lax hierarchies (13. 7p are completed by 



[d-Ai-qi{si + l),L^] ^a = 0- 



They differ from the rest of equations belonging to the dispersionful Whitham hierarchy as 
they are partially lattice systems in variables Sj. With supplementary equations fl3.19ll3.20ll 
the dispersionful Whitham hierarchy is in one to one corespondence to the original one. These 
equations are well defined as using the relation (13.161) we can pass from pseudo-differential op- 
erators to shift ones and formulate the dispersive analogue of the Whitham hierarchy entirely 
by means of shift operators. This is a counterpart of the alternative formulation performed 
in the Section 12.31 

We define the following Lie algebras of shift operators 



^a^AH© J^a^AH z = 1, 2, AT, 



1^0 ) K KO 



that are decomposed into Lie subalgebras with respect to the commutator. The different 
Lie algebras at are mutually connected through the relations fl3.17p and (13.181) due to which 
also the rules of expansion can be formulated. Thus, the algebras of pseudo-differential 
and shift operators act compatibly on the eigenfunctions i.e. Q{i,±)'^a = ci(i,±)^a and 
0(oo,±)^a = ci(i,±)^a5 for arbitrary i. Hence, the rules of expansion and i?-matrices related 
to decomposition of Lie algebras are invariant under the passage between both alternative 
formulations. Moreover, the form of Lax hierarchies (13.71) and zero-curvature equations (13.81) 
is preserved. Then, the set of linear equations (13. 5p following from (13. 2113. 3p is given by 




1=0 

and the definition of generating operators (13.41) remain unchanged. Thus, equation (13.161) 
follows from (13.60 for (3 = oo,n = 1. The zero-curvature equations (13. 8p together with 
Lax hierarchies (13. 7p defined for Lax operators (13.211) can be extracted from the 'charged' 
mult i- component KP hierarchy [13]. As we are using simultanously shift operators Aj with 
respect to different Si we call such a formulation of the dispersionful Whitham hierarchy as 
the mult i- component Toda-like. 

In the above formulation there is still missing the 'quantum' counterpart of Qko = In 
from Section 12. 3[ The natural logarithm of the shift operator can be defined in the way 
In Afc := -£;^k\a-o' ^^^^^ arbitrary. Over the logarithm In A^,, for fixed k, the algebra 
can be extended in such a way that the Lie algebra splitting is be preserved. Thus, we have 



d 

" fifco^a = InAfc^'^ = — e 



«=0 



Hence tko can be identified with Sk and the related equations from the dispersionful Whitham 
hierarchy, in the multi-component Toda-like formulation, are trivially satisfied as it should be. 
In the case of Lax hierarchies one obtains the translational symmetries (Lq,)^^^ = (Lq,)^^. As 
before, in the case of the original Whitham hierarchy, considering the bi-Hamiltonian formu- 
lation one would see that these translational symmetries generates higher order symmetries 
that are not included in our construction. The way of construction of the whole hierarchy 
generated by logarithmic operators, defined by means of dressing operators, for the two-field 
lattice Toda system is presented in [12j. 



Notice that In Afc^^ 7^ \n{d- qk{sk + 1)) as A^^a ^ {d - qk{sk + 1))° ^ a- Besides, 
the action of In {d — q{sk + 1)) on eigenfunctions would not be well defined. Nevertheless, 
one would like to have the formulation of the above logarithms by means of pseudo-differential 
operators. To calculate it one should use a formula similar to f l3.25p . but the calculations are 
not straightforward anymore. 

Example 3.4 To have a full dispersive counterpart of Example (12.1 p we have to complete the 
set of equations from Example (13. ip and calculate systems with the evolution parameters Sk- 
We will use the related operators (13. 4p from auxiliary linear equations (13. 6p in the form 

(^,)^ = (Afc + qk{sk + 1)) ^« (^a)^^ = (-fliAr') ^« 
(^„), = (A^ + (gfc(sfc + 2) + qk{sk + 1)) Ak + {qk{sk + 1))' + a{sk + 1) + a) 
Then, for t < — > yi and t < — > Si one finds equations (13. 9p and 

{qi)x = w - w{si - 1). 
For Si < — > Sj and yi < — > yj one finds systems 

qj - qj{si -l) = qi- qi{sj - 1) = (In \qi - qj\)^ (3.22) 

\ - (n\ _ ajajisj - 1) - aijsj - l)aj 
g,(.,-l)-5,(..-l) ' 

respectively. Finally, for Si < — > yj we have equations 

{qi)y, = aj - aj{si - 1) {ai)^ = ai{qi{si + 1) - qi) 

from which eliminating q.-^ one finds the multi- dimensional Toda equation. 

3.4 Finite-field reductions 

Now, in the representation of shift operators we are able to construct finite field reductions 
(13.1 op for arbitrary tIq, > given by a Lax operator C in the form 

(jioo— 2 N rii \ 

A^- + uA^--' + J2 ^i^i + E E ]^ = E^ (3.23) 

1=0 i=l k=l J 

being a dispersive counterpart of (I2.26p . where now u = Yl]=i Iki^k + j). Notice that some 
terms of (13. 4p can have coefficients that are nonlocal on lattice, see Example 13.61 However, as 
we are interested in the construction of (H-l)-dimensional equations we have to show how to 
represent (I3.23P only with A^ for fixed k. The equation (I3.22p belongs to the zero-cutvature 
equations that must be allways satisfied. Thus, from (I3.22p we have the relation 

qk{sk + l,Si - 1) - qi = qk{sk + 1) - qi{sk + !)• 

Hence, by (I3.18p one finds that 

A-^m = {Ak + qk{sk + 1, Si - 1) - qi)-' ^ = (Afe + qk{sk + 1) - qi{sk + 1))"' 

which when use recursively leads to 

Ar'"^ = (Afc + qkisi + 1 - m, + 1) - qi{si + l-m,Sk + 1))'^ ■ ... 

... ■ (Afc + qk{sk + 1) - qi{sk + l))~' ^ 

for m > 0. 



J- 



Example 3.5 The multi- component Toda hierarchy. 

Let us use the notation u{sk + m) := u'^'^\ The linear equation (13 .14^ represented by Afc 
transforms to 



C"^ = I Afc + qi^^ + ttkAk ^ + ^ai(^Ak + qi 



Then for ^ta^i = (^fc + ql )^ one calculates the multi- component Toda system 

(tti)^ = af ^ (^gf ^ - ^) + (^q^^^ - q,^ 
{<li)x = [qi - (qi - qk) + ^ - aj-" 



z.e. t/ie lattice analogue of (12.271) . For Vtki^ = — afcA^^\l/ we have the counterpart of (13.151) 
in the form 

,(-1) 



^ ^ _ (-1) 

i^i)yk ~ ^k 



/ N It qi - qk 



(1) (1) / ' 
Ti -Tk ' 



where i ^ k. 

Example 3.6 We will use the notation q := qi, A := Ai, y := yi and u^"^^ := u{si + m). 
The linear equation being a dispersive analogue of (12.281) is 

= {A + g(^) + uA-^^ + vA-^) ^. (3.24) 

Then, for fiooi* = (A + g^^)) ^ one finds 

Ux = u {(f^^ — g) + v^^^ — V 
V. = V (g(i) - g(-i)) 

gx 

analogue of (I2.29p . Aexi system is 



qx = u — 



Uy = a (g^^^ - g) 



= -aA ^ = na^ - ^^a 



qy = a — 



where a satisfies aa'^ = v. There are two simplest possible solutions on the field a: 

^(1)^(3)^(5) . ^^(-2)^(-4) . 

a = — 77:^r-FTs or a 



W(2)^;{4) . ... ^;{-l)t;{-3) . 



Example 3.7 /n order to complete the dispersionful version of Examples or \2.8\ we need 
to write the Lax operator C from Example \3.3\ in terms of the shift operator A := Ai as 

= (A^ + {qW + q(2))a + + U + aA-i) tfj. 

Now, from fiooi^P = (A + q^^^^ \I/ we find the lattice analogue of fl2.3UI) 

(g(l)+g(2))^ = „(l)_^ 

+ = -a 
a^ = a (g(^) - q) . 

On the other hand, by introducing b := q^^^ + q^'^^ and c := (q^^^^ + u, so that = 
(A^ + 6A + c + aA~^) we have also the lattice analogue of fl2.3ip 

k = a(') - a 

fioo2^ = (A^ + + c) ^ ^ Q = aWf>-afo(-i) 

at = a(c — c''"^-'). 

As we have been able with the use of expression (13.171) to construct Lax operators C, being 
dispersive counterparts of the meromorphic functions with 'finite' poles of higher order, we 
would like to use f l3.17p in the opposite way and represent £ by means of pseudo-differential 
operators. Thus, one finds the useful formula 

A--^ = {d- qkisk -m + ■ ... ■ {d - qkisk - l))-\d - qk)-'^, (3.25) 

where m > 0. However, shifted fields will appear, which for a = oo can be understood as new 
independent dynamical fields. So, one is able to construct closed equations being finite-field 
reductions (13.101) of KP hierarchy, also for nj higher then 1, that are analogues of algebraic 
orbit reductions of Whitham hierarchy. The above dispersionful analogues in general will 
have more dynamical fields then its dispersionless limits. The reductions of this kind of KP 
hierarchy have been considered earlier in [HI [9] in the context of multi-matrix models. 

Example 3.8 The linear equation (I3.24p written by means of pseudo- differential operators 
takes the form 

c^^J = {d + u{d - qy^ + v{d - qT\d - qy^) ^, 

where q' := q^~^\ Then, for i7oo2^ = (d^ + 2ii) \E' one finds the following system 

Ut = + ^v^ + 2{uq)^ 
Vt = V2x + 2vq^ + 2{vq')^ 
qt = -q2x + + 2gg^ 



I't = -qL - 2?2a' + + 2q'q 



'J 

X 



The evolutions of fields q and q' = q{si — 1) in the quasi- classical limit will be identical, 
i.e. four-field equation (13.260 is well defined dispersive counterpart of three-field (12.190 . This 
system was constructed in 

In the forthcoming article [26] the theory presented in this section will be treated more 
broadly and with all necessary proofs. The aim of [26] will be a construction of (1+1)- 
dimensional field and lattice soliton systems being dispersive counterparts for a wider class 
of dispersionless equations considered in [22j. 



4 Comments 

In [13] Takasaki provided a scheme for deriving the Whitham hierarchy as a dispersionless 
hmit of the 'charged' multi-component KP hierarchy. The charges are introduced through 
extra set of discrete variables Sj. The scheme is based on the r-f unction bilinear formal- 
ism [571 [2H] where the total charge Sqo + Sj equal to zero. In this sense the 'charged' 
multi-component KP hierarchy can be understood as a mult i- component generalization of 
Toda hierarchy. As a result, the scalar auxiliary linear equations (13.5113.61) . defined by means 
of shift operators, leading to the Lax formalism (13. 7113. 8p can be extracted. In literature 
there are several different formulations of the multi-component KP hierarchy, for example the 
Grassmannian approach [22] or the 9- method [30]. So, it seems that a further discussion of 
the corresponding connections between them and the dispersionful Whitham hierarchy from 
the present work and [131 E] is needed. 

Besides, due to the previous comments on the bi-Hamiltonian structures of the original 
and dispersionful Whitham hierarchy we see that the existence of the symmetries leads to 
the conclusion that these hierarchies are incomplete. So, we can extend this incompleteness 
onto multi-component KP and Toda hierarchies. This problem seems to be worth of further 
investigation. 
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